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Abstract. In this paper we provide a thorough classification of Friedman-Lemaltre- Robertson- 
Walker (FLRW) cosmological models in terms of the strong or weak character of their 
singularities according to the usual definitions. The classification refers to a generalised Puiseux 
power expansion of the scale factor of the model around a singular event. 



1. Introduction 

In the last few years the number of discussed matter contents for the different ages of our 
universe has increased. Besides the original dust, radiation and cosmological constant terms, 
researchers have come to include other contributions to the right-hand side of Einstein equations, 
being quintessence and phantom energy some of the most usual. The consideration of these new 
terms has been motivated by the attempts of explaining the experimentally inferred accelerated 
expansion of the universe [ij . 

The inclusion of these new matter contents, which do not satisfy some of the classical energy 
conditions, has as a consequence the appearance of new types of singularities which did not 
come up in former models, such as big rip [2] and sudden singularities [3], and other non-singular 
features, such as bounces or extremality events. All of these have been comprised under the 
name "cosmological milestones" by Cattoen and Visser [3] together with the classical big bag 
and big crunch singularities. 

In that paper the authors perform a classification of FLRW cosmological models in terms of 
the first coefficients of a generalised Puiseux expansion of the scale factor in time coordinate 
around one of these cosmological milestones. These coefficients are used to determine violations 
of energy conditions and appearance of polynomial scalar curvature singularities and derivative 
curvature singularities. Generically the classification depends on just the first three exponents 
of the Puiseux expansion. 

What we would like to do now is to complete the classification by considering other definitions 
of singularities. Geodesic incompleteness is commonly accepted as an indicator of the existence 
of a singularity in a space-time [5] and may happen even in cases where there is no polynomial 
scalar curvature singularity. 

Furthermore, it has been shown that even in the cases where a causal geodesic is incomplete, 
this does not mean that finite objects are necessarily crushed on approaching the singularity. 
These are considered weak singularities. In our classification we take into account this fact with 
the most common definitions of strong singularities. 



This involves calculation of causal geodesies in the corresponding space-times. This topic is 
reviewed in section [2l In section [3] geodesic equations are solved for the power expansion of 
the scale factor and the differentiability of the geodesies is analysed. Finally, in section H] the 
strength of the singularities, if any, is discussed in relation to the values of the exponents of the 
power expansion. The conclusions are summarized in section [5l More details about this issue 
may be found in [6]. 



2. Geodesic equations 

Generally calculation of geodesies in a space-time is a cumbersome task, since it requires solving 
a system of four ordinary quasilinear differential equations. Geodesies are parametrised by their 
proper time, 

dr^ = —gijdx^dx\ (1) 

where gij are the components of the metric tensor of the space-time in the coordinate chart 
provided by {xq, xi, X2, X3}. Proper time is defined up to a change of scale and origin, f = ar+b, 
and hence it is also called affine parameter. We denote by a dot derivatives with respect to proper 
time. 

A geodesic is said to be complete if it can defined for all values of r. On the contrary, it is 
said to be incomplete in the past (future) if it can be extended just to a value tq instead of —00 
(00). 

Geodesies are defined as the curves F on the space-time for which the length functional. 



L[T]= ds, ds' = -dr' = gijdx' dx\ (2) 



so 



has a extremum. 

The corresponding Euler-Lagrange equations, 



+ T]j,x^x^ = 0, (3) 

may be written in terms of the Christoffel symbols for the metric tensor, 

1 

= 2^*' {9ij,k + 9ik,j - gjk,i} , (4) 

plus an additional equation, 

S = -gijx'x^, (5) 

which simply states that we are using proper time as parameter. The constant 6 takes value one 
for timelike, zero for lightlike and minus one for spacelike geodesies. 

We use spherical coordinates, {t, r, ^,0}, with the usual ranges and t is coordinate time. 
Geodesies are therefore described providing {t{T),r{T),6{T), (P{t)) . 

In the case of FLRW eosmologieal models, 

ds^ = -df + a\t) { f\r)dr^ + (dO^ + sin^ ed(t>'') } 

/'W = rZl^' fc = 0,±l, (6) 

the large number of isometrics allows us a quick integration of geodesic equations. 

Since the space-time is homogeneous and isotropic, geodesies are straight lines and hence we 
may restrict the discussion to a plane 6* = 7r/2, ^ = const, choosing as origin of coordinates one 
of the points of the geodesic. 



Furthermore, changing the radial coordinate to 

{arcsinh r k = —1 ( sinh^ R 

r k = , ds"^ = -dt^ + a^{t) I dR^ + R^ (dO 

arcsin r k = \ I sin^ R 



2 , -2 

+ sm 




it is easy to check that Or = dr/ f{r) is another generator of isometries and therefore 

P = a\t)f{r)r (8) 



is conserved along geodesies. 

We are left then with just one equation, 



and the quadrature 

P 



which may be integrated after solving the equation in t and therefore need not be considered 
here. 

At this point it is clear that all information about the geodesies is encoded in the scale factor 
a{t). 

We consider just future-pointing geodesies, t > 0. 

Without losing much generality, we assume that the expansion, 

a(t) = co|t-to|* + ci|t-tir^ +••■, (11) 

where the exponents ry^ are real and ordered, 

7?o < ??i < • • • 

is valid close to a cosmological milestone at to- The coefficient cq must be positive in order to 
have a positive scale factor. 

We consider just singularities in the past, t > tQ,va. order to avoid signs and absolute values. 
Since equations are time-reversal symmetric, no information is lost with this restriction. 

At lowest order, ryo, in the flat universe case, /c = 0, the model behaves like a perfect fluid of 
density p and pressure p with a linear equation of state, 

2 , , 

p = wp, w = -1 + (12) 

Also at lowest order, three different behaviours are possible for the scale factor at to: zero, 
finite and divergent: 

• r/o > 0: the scale factor vanishes at to and generically we have a big bang or big crunch 
singularity. 

• r/o = 0: the scale factor is finite at to- If o^(0 is analytical, the event at to is regular. 
Otherwise a sudden singularity comes up [3l [7] . 

• r/o < 0: the scale factor diverges at to and a big rip singularity comes up. 

Since completeness of just causal geodesies is required for the analysis of singularities, we 
focus only on lightlike and timelike geodesies. 



Table 1. Derivatives of lightlike geodesies at to- 
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3. Geodesic completeness of causal geodesies 

Lightlike geodesic equations are straightforwardly integrated, 



a{t)i = P^ f a{t')dt' = P{t - To). (13) 



At lowest order, 



a(t) ~co|t-to!* ^t-to+ <( ■ (14) 

Ce^^/'^o % = -1 

Since generically t behaves as a power 1/(1 + r/o) of proper time, different levels of regularity 
appear depending on the value of 7]q. 

It is worth mentioning that for ijq lower or equal than minus one, lightlike geodesies do not 
reach the cosmological milestone at to, since it would take them an infinite proper time to reach 
it. They therefore do not see the singularity. This limiting case, which corresponds to a model 
with w = —5/3 when the universe is flat, has been named superphantom and considered in [8j. 

Results on the differentiability of lightlike geodesies at rjQ are consigned in table [TJ 

For the values of r/o for which the class of differentiability increases by one, it is necessary to 
consider further terms of the Puiseux expansion. 

We notice that, as rjQ decreases, the class of differentiability increases. 

Timelike geodesies may be analysed similarly, though in this case geodesic equations cannot 
be solved analytically. 
At lowest order, 



t=yi + ^:^^l + ^(t-to)-2^o, (15) 

the geodesic is singular at t = to for r/o > 0, since t blows up. 

However, the derivative is well deflned for negative (big rip) r/o- Near to we may write 



p2 



l + §^(t-to)-2^°, (16) 



Table 2. Derivatives of timelike geodesies at to- 
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and similar expressions for higher derivatives, 

- (t-to)"^''°""+^- (17) 

Again, as it is shown in table [21 the class of differentiability of timelike geodesies increases 
as r/o decreases. However, every timelike geodesic reaches the cosmological milestone at to and 
there are no curves with all finite derivatives there. 

4. Strength of singularities 

In the previous section we have shown that qualitatively the strength of singularities at to 
decreases with ijq, since the class of differentiability increases. We proceed now to check this 
qualitative statement with the most usual definitions of strong singularities. 

Roughly speaking, a singularity is considered strong if tidal forces are capable of disrupting 
a finite object falling into it 

This concept has been developed in several ways. For Tipler [TU] the finite volume is spanned 
by three Jacobi fields that form an orthonormal basis with the velocity u of the geodesic. The 
singularity if strong if the volume tends to zero on approaching the singularity. 

Another definition is due to Krolak [llj . for which the singularity is strong if the derivative 
of the volume is negative close to the singularity. Obviously, if a singularity is strong according 
to Tipler's definition, it is strong according to Krolak's, but not conversely. 

These definitions are meant to be used for gravitational collapse and therefore do not consider 
the possibility of big rip singularities. But these may be included in the framework just reversing 
a sign. 

Both definitions have been written in an amenable form by Clarke and Krolak [12j in terms 
of integrals of Riemann components along the geodesies. 

In our case the situation is even much simpler since the space-time is conformally flat and 
the Weyl tensor vanishes. 

For instance, a lightlike geodesic meets a strong singularity, according to Tipler's definition, 
at proper time tq if and only if 

r dr' r dT"Riju'u^ (18) 
Jo Jo 

diverges as r tends to tq. 



Table 3. Degree of singularity of null geodesies around to- 
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And a lightlike geodesic meets a strong singularity at proper time tq if and only if 

dr'Riju'u^ (19) 



diverges as r tends to tq. 

In our case, the velocity of the geodesic is 



R v^ui - 2P2 ('^ - ^ ] ~ ''-^ + + . . . (20) 



and at lowest order there are two cases, depending on whether the curvature term dominates 
over the first term. 

The results according to both definitions are summarized in table [3l 

As we see, besides the models with rjQ < —1, for which no lightlike geodesic reaches the 
cosmological milestone at to, there are only two cases without strong singularities: some of 
the cases with rjo = 0, which are named sudden singularities, some of which had already been 
studied in [13]; and some of the cases with ijq = 1, cq = 1, k = —1, which are at first order Milne 
universe, which is Minkowski empty space in other coordinates. 

Again, there are limiting cases which require resorting to further terms in the expansion. 

As it was pointed out, it is explicitly checked that Krolak's definition includes more cases 
than Tipler's. 

The analysis of the strength of singularities of timelike geodesies is somewhat more involved, 
since there are no both necessary and sufficient conditions for the appearance of strong 
singularities: 

According to Tipler's definition, a timelike geodesic meets a strong singularity at proper time 
To if 

r dr' r dT"RijU^u^ (21) 
JO Jo 

diverges as r tends to tq. 

With Krolak's definition, a timelike geodesic meets a strong singularity at proper time tq if 

dr'Riju'u^ (22) 



Table 4. Degree of singularity of the fluid congruence of timelike geodesies around Iq. 
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diverges as r tends to tq. 

Necessary conditions are slightly different. With Tipler's definition [121, ^ timelike geodesic 
meets a strong singularity, then 

/j(r) = l^dr'l^^ dr"\Rl^iu'u'\, (23) 

diverges as r tends to tq for some i, j, where the components are referred to a parallely 
transported orthonormal frame. 

With Krolak's definition, if a timelike geodesic meets a strong singularity, then 

/Kt)= fdr' i?L,n^' , (24) 







diverges as r tends to tq for some i, j. 

Fortunately, this set of sufficient and necessary conditions is accurate enough to allow a 
thorough classification of the singularities of timelike geodesies in FLRW models. 

There are two sets of timelike geodesies with different behaviour: 

Timelike geodesies with P = 0, for which the time coordinate is essentially proper time, 

t-to = T-To, r = ro, (25) 

form the congruence of fluid worldlines, since the coordinates are comoving for the perfect fluid, 
u = dt, and do not suffer therefore any problems of differentiability. On applying the conditions 
for the appearance of strong singularities, we reach the results of table HI 

Along timelike geodesies with radial velocity, P ^ 0, strong singularities appear in more 
cases, as we show in the results comprised in table [5l 

As we see, the results are essentially the same as for lightlike geodesies, with a difference of 
behaviour for models with coeficient t]q lower or equal than minus one. Timelike geodesies do 
reach the cosmological milestone at to in the form of a strong singularity. These models are 
lightlike geodesically complete, though they are timelike geodesically incomplete. In fact, table 
provides the classification of FLRW cosmological models according to the strength of their 
singularities. 

5. Conclusions 

In this paper we have provided a complete classification of FLRW cosmological models according 
to the strength of their singularities in terms of a generalised Puiseux expansion of the scale 
factor in coordinate time around cosmological milestones. 



Table 5. Degree of singularity of timelike geodesies around io- 
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Though the velocity of the geodesic is finite at big rips, this does not prevent the appearance 
of strong singularities, except for two groups of models: those which behave at lowest order as 
Milne universe and those with sudden singularities. However, lightlike geodesies do not reach 
the big rip singularities for exponents rjQ lower or equal than minus one. 

The results of Cattoen and Visser in [3] are more restrictive, since they do not deal with 
the strength of the curvature singularities. For them models with % = 0, r/i > 2 or ryi = 1, 
r]2 > 2, and those with r/o = l,A; = — 1, co = l, ?/i>3 are free of polynomial scalar curvature 
singularities. These results coincide with ours for Milne-like models with Krolak's definition for 
strong singularities, but include more models in the case of sudden singularities. 

Their results for derivative curvature singularities are even more restrictive. The only models 
which are free of such singularities are those with ryo = and natural exponents iji, i > 1 and 
those with rjo = 1, k = —1, cq = 1 and natural exponents ?/i > 3, i > 1. 

This is not surprising, since derivatives of the curvature tensor do not appear in our equations. 

However, the main difference between both results arises from the fact that curvature 
singularities do not see that lightlike geodesies do not reach the singularity in finite proper 
time. 

Though this classification of singularities conveys the idea of ubiquous singularities in FLRW 
cosmological models, it is worthwhile mentioning that singularities mostly appear in models 
with vanishing, divergent or non-smooth scale factors. 

A similar scenario appeared in inhomogeneous scalar field Abelian diagonal G2 models |14j . 
where singularity-free cosmological models formed an open set. 
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